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Abstract 

We show that the median m{x) in the gamma distribution with 
parameter a; is a strictly convex function on the positive half-line. 
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1 Introduction 

The median of the gamma distribution with (positive) parameter x is defined 
implicitly by the formula 
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e -H x - 1 dt = =- / e-H x - l dt. (I) 
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In a recent paper (see (Sj) we showed the < m'(x) < 1 for all x > 0. 
Consequently, m{x) — x is a decreasing function, which for x = 1,2,... 
yields a positive answer to the Chen-Rubin conjecture. Other authors have 
solved this conjecture in its discrete setting (see [2], P, [3]). 

In [3] convexity of the sequence m(n + 1) has been established, and the 
natural question arises if m{x) is a convex function. The main result of this 
paper is the following. 

Theorem 1.1 The median m{x) defined in 0j satisfies m"{x) > 0. In 
particular it is a strictly convex function for x > 0. 
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2 Proofs 



The proof is based on some results in which we briefly describe. Con- 
vexity of m is studied through the function 

ip(x) = log—^-, x>0. (2) 

This function played a key role in [5], and we recall its crucial properties in 
the proposition below. 

Proposition 2.1 The function x — > ary(x) is strictly decreasing for x > 
and 

lim x(p(x) = log 2, 
x— >o+ 

lim xip(x) = —. 

x—+oo 3 

Remark 2.2 Provosition \2. 1\ is established by showing (xip(x))' < 0. It fol- 
lows that the function (p(x) is itself strictly decreasing and <p(x) < —X(p'(x). 

The starting point for proving Theorem II .11 is the relation [3 (10)] 

(x<p(x))' = -e 9(x \A(x) + B(x)), (3) 

where 

g{x) = x(<p{x) - 1 + e"^), 

e - v(1 -e-/-) (l - (l + £) e -/«) dSj 

= - / te-^C* + 1) 
^ Jo 

and where £ is a certain positive, increasing and concave function on [1, oo) 
satisfying £'(t + 1) < 8/135 for t > 0. To establish these properties of £ is 
quite involved, and we refer to Section 5] for details. 

Before proving the theorem we state the following lemmas, whose proofs 
are given later. 

Lemma 2.3 For the function g we have 

g(x) < xip(x), 

—g{x) < — x<p (x)ip(x) and 

—g'{x) < —xip'(x) 

for all x > 0. 
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Lemma 2.4 We have for x > 

x(p(x) 3 



A(x) < 
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Lemma 2.5 W^e have for x > 

4 



Proof of Theorem M.ll From equation (J2J) we get 

m»(x) = -e-vW (2tf/(x) + x<p"{x) - xtp'(x) 2 ) , 
so that m"(x) > is equivalent to the inequality 

(x(p(x))" < xip'(x) 2 . 

Differentiation of © yields 

(x<p(x))" = e^ x \-g'{x)){A{x) + B{x)) + e 9 ^(-A'{x) - B'{x)). 

By using Lemma 12 . 41 and 12 . 51 it follows that 

1 8 1 

-A'{x)-B\x) < --x<ff{x)<p(x) 2 + ^3 - -<p{xf 



1 I, \ f si <f( X ) ( 8 1 / » 
■—XCp (Xjip(x) H ; : — 7TTT XifiiX 

2 Yy) ^ y > x \l35(x<p(x)) 2 6 ^ V ' 



Here the expression in the brackets is positive, since (xip(x)) s < (log2) 3 < 
48/135. Therefore, and because <p(x) < —x(p'(x), 



-A'[x) - B'(x) < ^xip'(x) 2 xip(x) + xif'(x) 2 f 



135(x(^(x)) 2 6 



-xtp(x) 



x ^ {x)2 [ mxl(xw + 1 3 Mx: 



We also have from Lemma 12.31 12.41 and 12.51 

-g'(x)(A(x)+B(x)) < + 



< X (f (x) (p{x) I — h 



6 135(x<f(x)y 
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Combination of these inequalities yields 
Supposing that x > 1, it follows that 

= x(^'(x) 2 /ii(x(^(x)), 
where h\ is given by 

One can show that /ii attains its maximum on the interval [1/3, log 2] at 
the left end point and that /ii(l/3) = (551/810) « 0.9494. Therefore it 
follows that (x(p(x))" < x(p'(x) 2 for x > 1. 

For < x < 1 the estimate —g'(x) < —x<p'(x) from Lemma 12.31 is used 
and in this way we get 

(x<p(x))" < x<f'(x) 2 h2(x(p(x)), 

where 

Wi) = e * — h 1 + - . 

Since x < 1 and xip{x) decreases we have xcp(x) > <p(l) = — log log 2. One 
can show that attains its maximum on the interval [— log log 2, log 2] for 
t = —log log 2 and that hi{— log log 2) ~ 0.9616. Therefore (x<p(x))" < 
xip'(x) 2 for x < 1. □ 

Remark 2.6 The function hi becomes larger than 1 on the interval [1 /3, log 2] , 
so h>2 cannot be used to obtain the inequality (xip(x))" < xp'(x) 2 for all 
x > 0. 
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Proof of Lemma \2. 31 It is clear that g(x) < x<p(x). Differentiation yields 

-g'(x) = - V (x) + (l-x<p'(x))(l-e-^ x) ) 

< —y>(x) + (1 — x<p' (x))ip(x) = —xip' {x)ip{x) , 

where we have used 1 — e~ a < a for a > 0. 

To find an estimate that is more accurate for x near we use 

-g'(x) = -x(f/(x)(l-e-vto)-(p(x) + l-e-rt x ') 
< -x(p'(x) - <p(x) + 1 - e~^ x) < -x<p'(x). 

□ 

Proof of Lemma \2~^\ Using that 1 — e~ a < a and 1 — (1 + a)e~ a < a 2 /2 
for a > 0, we can estimate A{x) by 



rx<fi(x) / 2 \ 
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A computation shows that 



-A'(x) = -(if(x) + xy'{x))e- x ^ x) e x{l - e ~* ix)) {l - (1 + ip(x)) e -^) 
-j MX) e-^-^ (l-(l + g e -^) 2 d S 

rx<p(x) . „2 

+ / e s e ^-e- a/X ) S e s/* ds 

Jo x 

1 /-x(p(x) ^ 

< -{(p(x) +xip'(x))-ip(x) 2 + s 2 e~ s/x ds^r 

2 Jo x 6 

< -\^{xf -\x^'{x)^{x) 2 + ^f{xf 
= -^f( x f - \xip' {x)if{x) 2 . 

□ 

Proof of Lemma \2. 51 These estimates follow directly from the inequality 
£'(* + 1) < 8/135. □ 
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